2, 2008 1:21 WSPC/INSTRUCTION FILE jrpmplahepph 



Modern Physics Letters A 

© World Scientific Publishing Company 



Light scalars as tetraquarks or two-meson states from large N c 
and unitarized Chiral Perturbation Theory. 

J.R. PELAEZ 

Departamento de Fisica Teorica II. Universidad Complutense. 28040 Madrid. SPAIN 

jrpelaez@fis. ucm. es 

Received (Day Month Year) 
Revised (Day Month Year) 

By means of unitarized Chiral Perturbation Theory it is possible to obtain a remarkable 
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1. Introduction 

We review here how the spectroscopic nature of the lightest scalar resonances can be 
obtained from the large N c behaviofl] of their associated poleJ^in the meson- meson 
scattering one-loop Chiral Perturbation Theory unitarize cP amplitudes. 

The existence and nature of the lightest scalar resonances is a longstanding con- 
troversial issue for hadron spectroscopy. Their relevance is twofold: First, one of the 
most interesting features of QCD is its non-abelian nature, so that gluons interact 
among themselves and could produce glueballs, which are isoscalars. Moreover, the 
lightest ones are expected to be also scalars. Naively, once all quark multiplets are 
identified in the scalar-isoscalar sector, whatever remains are good candidates for 
glueballs. Unfortunately, the whole picture is more messy due to mixing, since the 
resonances we actually see are a superposition of different states. Second, it is also 
understood that QCD has an spontaneous breaking of the chiral symmetry since 
its vacuum is not invariant under chiral transformations. The study of the scalar- 
isoscalar sector is relevant to understand the QCD vacuum, which has precisely 
those quantum numbers. 
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Although QCD is firmly established as the fundamental theory of strong inter- 
actions and its predictions have been thoroughly tested to great accuracy in the 
perturbative regime (above 1-2 GeV), it becomes non-perturbative at low energies, 
and helps little to address the existence and nature of light scalar mesons. Hence, 
the usual perturbative expansion in terms of quarks and gluons has to be abandoned 
in favor of somewhat less systematic approaches in terms of mesons. As a matter of 
fact, most chiral descriptions of meson dynamics do not include quarks and gluons 
and are hard to relate to QCD, and the spectroscopic nature is thus imposed from 
the start. In contrast, models with quarks and gluons, even those inspired in QCD, 
have problems with chiral symmetry, small meson masses, etc... 

An exception is the formalism of Chiral Perturbation TheorjEEl (Ch.PT), which 
is the most general effective Lagrangian made out of pions, kaons and etas, that 
respects the QCD chiral symmetry breaking pattern. These particles are the QCD 
low energy degrees of freedom since they are Goldstone bosons of the QCD sponta- 
neous chiral symmetry breaking. For meson-meson scattering ChPT is an expansion 
in even powers of momenta, 0(p 2 ), 0(p 4 )..., over a scale A x ~ 47r/o ~ 1 GeV. Since 
the u, d and s quark masses are so small compared with A x they are introduced 
as perturbations, giving rise to the n,K and r\ masses, counted as 0(p 2 ). At each 
order, ChPT is the sum of all terms compatible with the symmetries, multiplied 
by "chiral" parameters, that absorb loop divergences order by order, yielding finite 
results. ChPT is thus the quantum effective field theory of QCD, and it allows for 
a systematic and model independent analysis of low energy mesonic processes. 

Since ChPT is an expansion in momenta and masses, it is limited to low energies. 
As the energy grows, the ChPT truncated series violate unitarity. Nevertheless, in 
recent years ChPT has been extended to higher energies by means of unitarization 
HWrtmilal The results are remarkab Jf extending the one-loop ChPT calcula- 
tions of two body n, K or n scattering up to 1.2 GeV, but keeping simultaneously 
the correct low energy expansion and with chiral parameters compatible with stan- 
dard ChPT. Furthermore it generates the poles associated to the p(770), if* (892) 
and the octet 4> vector resonances, as well as those of the controversial scalar states, 
namely: the a, k, ao(980), and the /o(980). 

Among these states, the most controversial are the a, now called /o(600) in the 
Particle Data Group (PDG) Review El, and the n. They appear as broad resonant 
structures in the scalar channels in meson- mes on scatte ring since they do not display 
a Brcit-Wigner shape. Still, 

many group J7l*H 1H«|1 h| 

do find an associated pole in 
the amplitude, but deep in the complex plane. For a compilation of a and k poles 
see the nice overview in 1^1 Let us remark that meson-meson scattering datd^ 
are hard to obtain, since they are extracted from reactions like meson-V — >meson- 
meson-V, but with assumptions like a factorization of the four meson amplitude, 
or that only one meson is exchanged and that it is more or less on shell, etc... 
All these approximations introduce large systematic errors. Very recently, however, 
other experiments on meson-meson interactions have become available, like t he v ery 
precise determination of a combination of hit phase shifts from Kn decays ^1 or 
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the results from charm decays The latter have analyzed the pole structure of 
their amplitudes and seem to find both the a and n poles in reasonable agreement 
with the works mentioned above but in completely different processes. 

The controversy about their spectroscopic nature is even stronger. The interest 
of unitarized ChPT is that these states are generated from chiral symmetry and 
unitarity, without any prejudice toward their existence or structure. Amazingly the 
nine scalars are generated together, suggesting they form an SU(3) multiplct, with 
a similar composition, although probably mixed with other states when possible. 

Once again, ChPT has the advantage that it is the QCD low energy effective 
theory, and should beh ave as such in certain limits. For this reason we have studied^ 
the large N c expansiorPU, which is the only analytic approximation to QCD in the 
whole energy region. Remarkably, it provides a clear definition of qq states, that 
become bound states when jV c — > oo, whereas tetraquark or two meson states do 
not. Starting from the unitarized ChPT meson-meson scattering, and scaling the 
parameters with jV c according to the QCD rules^, it has been possible to show that 
vector mesons follow remarkably well their qq expected behavior, whereas all scalar 
nonet candidates behave as if their main component had a tetraquark or two-meson 
structure. 

In the next two sections we review in more detail the description of meson- 
meson scattering with unitarized ChPT. In Section 4 we determine the parameters 
of the poles in those amplitudes, and in Section 5 we will determine how those poles 
behave in the large N c limit. Finally we will present some brief conclusions. 



2. Meson-meson scattering and Chiral Perturbation Theory 

Chiral Perturbation Theory*^ (ChPT) is built as the most general derivative ex- 
pansion of a Lagrangian containing only pions, kaons and the eta. These particles 
are the Goldstone bosons associated to the spontaneous chiral symmetry breaking 
of massless QCD. In practice, and for meson-meson scattering, ChPT becomes an 
expansion in even powers of energy or momenta, denoted as 0(p 2 ),0(p 4 ), etc., 
over the scale of the spontaneous breaking, i.e.,4-7r/o ~ 1.2 GeV. Of course, quarks 
are not massless, but their masses are small compared with typical hadronic scales, 
and they are introduced as perturbations in the same power counting, giving rise to 
the 7r, K and rj masses, counted as M = 0(p 2 ). The main advantage of ChPT is that 
it provides a Lagrangian that allows for true Quantum Field Theory calculations, 
and a chiral power counting to organize systematically the size of the corrections 
at low energies. For example, in the isospin limit, the leading order Lagrangian is 
universal since it only depends on /o, the pion decay constant in the chiral limit, 
and the leading order masses Mj, Mjjr and M2. However, it is possible to calculate 
meson loops, whose divergences are renormalized in a finite set of chiral parameters 
at each order in the expansion. The dependence on the QCD dynamics only comes 
through higher order parameters. For instance, meson-meson scattering to one loop 
depends on eight Li parameters^, whose values are shown in Table 1. As usual after 
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renormalization, they depend on a regularization scale /z: 

LiifMt) = Li(m) + lo s— > (!) 

where F, are constant^. Of course, in physical observables the „ dependence is 
canceled through the regularization of the loop integrals. This procedure can be 
repeated obtaining finite results at each order. As long as we remain at low energies, 
only a few orders are necessary and the theory is predictive, since once the set of 
parameters up to that order is fixed from some experiments, it should describe, to 
that order, any other physical process involving mesons. 



Table 1. 0(p 4 ) chiral parameters (xlO 3 ) and their N c scaling. In the ChPT column, 
Li, Z/2i L3 come from 23 and the rest from 4 . The IAM columns correspond to different fits 3 



o( P 4 ) 


N c 


ChPT 


IAM I 


IAM II 


IAM III 


Parameter 


Scaling 


H = 770 MeV 


H = 770 MeV 


H = 770 MeV 


ft = 770 MeV 


Li 


0(N C ) 


0.4 ± 0.3 


0.56 ±0.10 


0.59 ±0.08 


0.60 ± 0.09 


L 2 


0{N C ) 


1.35 ±0.3 


1.21 ± 0.10 


1.18 ± 0.10 


1.22 ±0.08 


L 3 


0(N C ) 


-3.5± 1.1 


-2.79 ±0.14 


-2.93 ±0.10 


-3.02 ± 0.06 


L4 


O(l) 


-0.3 ±0.5 


-0.36 ±0.17 


0.2 ± 0.004 


(fixed) 


Ls 


O(iVc) 


1.4 ±0.5 


1.4 ± 0.5 


1.8 ±0.08 


1.9 ± 0.03 


L e 


O(l) 


-0.2 ±0.3 


0.07 ± 0.08 


0±0.5 


-0.07 ± 0.20 


L 7 


O(l) 


-0.4 ± 0.2 


-0.44 ±0.15 


-0.12 ±0.16 


-0.25 ± 0.18 


L 8 


0(N C ) 


0.9 ±0.3 


0.78 ±0.18 


0.78 ± 0.7 


0.84 ±0.23 


2Li - L 2 


O(l) 


-0.55 ±0.7 


0.09 ±0.10 


0.0 ±0.1 


-0.02 ± 0.10 



Other salient features of ChPT are its model independence and that it has been 
possible to obtairpE2 from QCD the Li large N c behavior, given in Table 1. Also, 
these parameters contain information about heavieJ^ meson states that have not 
been included as degrees of freedom in ChPT. 

However, ChPT is limited to low energies, since the number of terms allowed 
by symmetry increases dramatically at each order, because the ChPT series violate 
unitarity as the energy increases, and finally, due to resonances that appear rather 
soon in meson physics. These states are associated to poles in the second Riemann 
sheet of the amplitudes that cannot be accommodated by the series of ChPT, which 
are polynomial (there are also logarithmic terms, irrelevant for this discussion). 

For the above reasons, in rece nt years ChPT it has been extended to higher 
energies by means of unitarizatkm EE EEni]in B| that we discuss next. 

3. Unitarized Chiral Perturbation Theory 

In order to compare with experiment it is customary to use partial waves tjj of 
definite isospin / and angular momentum J. For simplicity we will omit the /, J 
subindiccs in what follows, so that the chiral expansion becomes t ~ t% + t± + .... 
with t% and £4 of 0(p 2 ) and 0(p 4 ), respectively. The unitarity relation for the partial 
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waves t{j, where i,j denote the different available states, is very simple: when two 
states, say " 1" and "2" , are accessible, it becomes 

ImT = TET* Imr _1 = -E => T = (ReT -1 - i E) -1 (2) 



with 



where <jj = 2qi/y/s and qi is the CM. momentum of the state i. The generalization 
to more than two accessible states is straightforward in this notation. Note that, 
since ImT -1 is fixed by unitarity we only need to know the real part of the Inverse 
Amplitude. Note that Eq.(J2J is non- linear and cannot be satisfied exactly with a 
perturbative expansion like that of ChPT, although it holds perturbatively, i.e, 

Im T 2 = 0, Im T A = T 2 E T* + 0{p 6 ). (4) 

The use of the ChPT expansion ReT -1 ~ T 2 _1 (l - (ReT 4 )T 2 _1 + ...) in eq.Q, 
guarantees that we reobtain the ChPT low energy expansion and that we are taking 
into account all the information included in the chiral Lagrangians (both about N c 
and about heavier resonances). In practice, all the powers of l//o in the amplitudes 
are rewritten in terms of physical constants f v or fx or f„. At leading order this 
difference is irrelevant, but at one loop, we have three possible choices for each 
power of /o in the amplitudes, all equivalent up to 0(p 6 ). It is however possible to 
substitute the /o's by their expressions in terms of / w or fx or f n in such a way 
that they cancel the 0(p 6 ) and higher order contributions in eq.QJ, so that 

Im T 2 = 0, Im Ti = T 2 E T 2 * . (5) 

We will call these conditions "exact perturbative unitarity". From eqs.J2Jl,©i an d 
the ReT" 1 ChPT expansion, we find 

T ~ T 2 (T 2 — T 4 ) -1 T2, (6) 

which is the coupled channel IAM that has been used to unitarize simultaneously 
all the one-loop ChPT meson-meson scattering amplitude^. Since all the pieces 
are analytic it is straightforward to continue analytically to the complex s plane 
and look for poles associated to resonances. Indeed, the analytic continuation to 
the complex s plane has also been justified in terms of dispersion relations in the 
elastic cas Alternative methods have been pro pose d and applied successfully to 
the one loop ChPT, for instance for irn scattering^. However, in this brief review 
we concentrate just on the IAM due to its remarkable success and simplicity, since 
it only involves algebraic manipulations on the ChPT series. 

The IAM was applied first for partial waves in the elastic region, where a single 
channel is enough to describe the data. This approach was able to generate the 
p and a poles in tttt scattering and that of K* in irK — * ttK. Later on, it has 
been noticed that the k pole is also obtained within the single channel formalism. 
Concerning coupled channels, since only the the tttt, ttK and nrj amplitude J^^ l were 
known at that time, at firslP' only the leading order and the dominant s-channel 
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loops were considered, neglecting crossed and tadpole loop diagrams. Despite these 
approximations, a remarkable description of meson-meson scattering was achieved 
up to 1.2 GeV, generating the poles associated to the p, K*, / , a , a and n. The 
price to pay was, first, that only the leading order of the expansion was recovered 
at low energies. Second, apart from the fact that loop divergences were regularized 
with a cutoff, thus introducing an spurious parameter, they were not completely 
renormalized, due to the missing diagrams. Therefore, it was not possible to compare 
the Li, which encode the underlying QCD dynamics, with those already present in 
the literature. In addition, the study of the large N c limit is cumbersome due to the 
incomplete renormalization and since we do not know how the cutoff should scale. 

As already explained the whole approach is relevant to study the existence and 
properties of light scalars and it is then very important to check that these poles 
and their features are not just artifacts of the approximations, estimate the uncer- 
tainties in their parameters, and check their compatibility with other data. These 
considerations triggered the interest in calculating and unitarizing the remaining 
meson-meson amplitudes within one-loop ChPT. Hence, the KK — > KK calcula- 
tion was completed^, thus allowing for the unitarization of the tttt, KK coupled 
system. There was a good agreement of the IAM description with the existing Lj. 
reproducing again the resonances in that system. More recently, we have completed 
the one-loop meson-meson scattering calculation, including three new amplitudes: 
Krj — ► Krj, rjrj — ► 7777 and Kir — ► Kr], but recalculating the other five amplitudes, 
unifying the notation, ensuring " exact perturbative unitarity" , Eq.j5). and also cor- 
recting some errors in the literature. Next, we have unitarized these amplitudes with 
the coupled channel IAM, which allows for a direct comparison with the standard 
low-energy chiral parameters, in very good agreement with previous determinations. 
In that work we presented the full calculation of all the one-loop amplitudes in di- 
mensional regularization, and a simultaneous description of the low energy and the 
resonance regions. 

The first check was to use the standard ChPT parameters, given in Table 1 
to see if the resonant features were still there, at least qualitatively, and they are. 
Thus, they are not just an artifact of the approximations and of the values chosen 
for the parameters. As already commented, this comparison can only be performed 
now since we have all the amplitudes renormalized in the standard MS — 1 scheme. 

After checking that, we made an IAM fit. Systematic uncertainties are the largest 
contribution to the resulting error bands as well as in the fit parameters in Table 
1. These error bands are calculated from a MonteCarlo Gaussian sampling^ (1000 
points) of the L.; sets within their uncertainties, assuming they are uncorrelated. 
Note that in Table 1 we list three sets of parameters for the IAM fit, fairly compatible 
among them and with those of standard ChPT. They correspond to different choices 
when reexpressing the fo parameter of the Lagrangian in terms of physical decay 
constants. The IAM I fit was obtained*" using just f„, which is simpler but unnatural 
when dealing with kaons or etas. There, it could be observed that the /o(980) region 
was not very well described yielding a too small width for the resonance. 
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Fig. 1. I AM fit to meson-meson scattering data, set II in Table 1. The uncertainty bands are 
dominated by systematic errors in the data 17 . 



For that reason, Fig.l shows the results of a second fit (IAM II) using ampli- 
tudes written in terms of fx and f v when dealing with processes involving kaons 
or etas. Let us remark that the data in the /o(980) region is well within the uncer- 
tainties. In particular, in the 0(p 2 ) one factor of l// ff has been replaced^ by 1/ fx 
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for each two kaons present between the initial or final state, or by 1//^ for each two 
etas appearing between the initial and final states. In the special case Kn Kit 
l/f% has been changed by l/ifxfr])- The difference between the two ways of writ- 
ing the leading order amplitudes is 0(p 4 ), a nd is therefore included in the next to 
leading order contribution using the relation^ between the decay constants and 
fo. The l//o factors in each loop function at 0(p 4 ) (generically, the J(s) given in 
the appendix of™) are changed to satisfy "exact perturbative unitarity", eqs.|(SJ- 
According to the ChPT counting, the amplitudes are the same up to 0(p 4 ), but 
numerically they are slightly different. A similar choice has been suggested™^ inde- 
pendently within non-unitarized standard SU(3) ChPT to avoid the uncertainties 
arising from fluctuation of vacuum ss pairs. In particular, it is suggested to obtain 
elastic amplitudes Apq from the safe combinations F^A^ or F^F^A^k, etc... For 
external fields this amounts to our choice, and the normalization of internal loop 
function is then dictated by exact perturbative coupled channel unitarity. From Ta- 
ble 1, we see that the only sizable change is in the parameters related to the decay 
constants, i.e., L4 and L5. For illustration we give in Table 1 a third fit, IAM III, 
obtained as IAM II but fixing L4 = 0, its large N c limit, as it is done in recent K14 
0(p 4 ) determinations. As seen in Table 1, these chiral parameters are compatible 
with those from standard ChPT, which means that we have a simultaneous descrip- 
tion of the low energy and resonance regime. Finally, since the expressions are fully 
renormalized, all the QCD N c dependence appears correctly through the Li and 
cannot hide in any spurious parameter. 

At this point one might wonder about higher order effects. There is a simple way 
to extend the IAM to higher orderJ^, first applied to two loop irir scattering and the 
pion form factoJSi, with remarkable results. This study has been extended^ with 
a careful analysis of the uncertainties. The amplitude depends on the 0(p 4 ) and 
0(p 6 ) parameters through six combinations, called &j. Despite the poor knowledge 
about these two-loop parameters a good description of the data is founds, including 
the a and p regions, with p aram eters compatible within errors with those in the 
literature. The error analyst is also of relevance because it was shown that the 
IAM crossing violations, taking into account the present experimental uncertainties, 
are " not very large in percentage terms" . 



4. Poles associated to resonances 

The most interesting feature of chiral unitary approaches is that the poles thus 
generated are not included in the original ChPT Lagrangian and hence appear 
without theoretical prejudices toward their existence, classification in multiplets, or 
nature. Remarkably, the scalar resonances a, k, ao(980) and Jo (980) appear together 
in those chiral unitarized amplitudes, and it seems natural to interpret them as a 
nonet (see alscEHI). Nevertheless, we should dis ting uish two different resonance gen- 
eration mechanisms: it had already been notecM™, that to generate the scalars just 
the leading order and a cutoff was enough, whereas the vector mesons require the 
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chiral parameters, particularly^l&l L\, L2 and L3. Of course, the chiral parameters 
are always present, but changes within errors in their values do not affect the ex- 
istence of the light scalar poles, but just the details of their description. Since the 
vectors are fairly well established qq states, this difference suggests that scalars like 
the cr, k, etc, may have a different nature. With the amplitudes described in the 
previous sections we expect to reach a more conclusive statement. 

Thus, in Table 2 we show the pole position for the resonances^, including uncer- 
tainties, for the different IAM parameter sets given in Table 1. For comparison we 
also provide those obtained in the "approximated" whereas we list in Table 

3 the poles listed presently in the PDcPl. These results deserve some comments: 

• The vectors p(770) and if* (892), are very stable within chiral unitary approaches. 
Their positions are almost the same irrespective of whether one uses the single 
channel^, the approximated coupled channepl or the complete IA]\P. 

• The a and k pole positions are robust within these approaches. No matter what 
version of the IAM is used. Note the small uncertainties in some of their param- 
eters, in very good agreement with recent experiments^! 

• The /o(980) has a sizable decay to two different channels and therefore it can only 
be studied in a coupled channel formalism. In practice, both the approximated 
and complete IAM generate a pole associated to this state at approximately 
the same mass. However, as already remarked, the unitarized ChPT amplitudes 
using juslPl /jr-yield a too narrow width. The good news is that it can be well 
accommodated^ using /„■, fx and /,,. 

• The ao (980) also requires a coupled channel formalism, and the data on this region 
is well described either by the approximated or the complete IAM. However, 
the presence of a pole is strongly dependent on whether we write the ChPT 
amplitudes only in ter ms o f f n , or using the three f n , fx and f n . In the first case, 
it has been pointed ou that the use of the "approximated" IAM with just /„., 
favored a "cusp" interpretation of the ao(980) enhancement in irr) production. 
With the complete IAM we do not find a pole near the ao(980) enhancement 
and indeed the nr/ phase-shift does not cross 7r/2 and it neither has a fast phase 
movement. In contrast, when expressing /o in terms of f n , ]"k and f n as described 
in the previous section, we do find a pole and its associated fast phase movement 
through 7r/2, either with the approximated or complete IAM. Thus, this pole is 
rather unstable as can be noticed from its large uncertainties in Table 2. 

• We are more favorable toward the pole interpretation for the do (980) because, 
first, it is also able to describe better the /o(980) width. Second, we have already 
remarked that even within standard SU(3) ChPT, the scattering amplitude calcu- 
lations are safer against ss vac uum fluctuations if these fluctuation terms, related 
to L4 and Lq are absorbed^-*- 1 in terms of physical decay constants , fx and /, ; 
when expanding the anplitudes. 

Let us remark that the /o(980) and ao(980) resonances are very close to the KK 
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Table 2. Pole positions (with errors) in meson-meson scattering. When close to the real axis 
the mass and width of the associated resonance is y/s po i e ~ M — iF/2. 



-/Sp o!e (MeV) 


P 


K* 


a 


fo 


ao 


K 


1AM Approx 
(no errors) 


759-i 71 


892-i21 


442-i 227 


994-i 14 


1055-i21 


770-i 250 


IAM I 


760-i82 


886-i21 


443-i217 


988-i 4 


cusp? 


750-i226 


(errors) 


± 52± i25 


± 50± i8 


± 17± 112 


± 19± i3 




±18±ill 


IAM II 


754-i 74 


889-i 24 


440-i212 


973-ill 


1117-i 12 


753-i 235 


(errors) 


± 18± ilO 


± 13± i4 


± 8± i 15 


+39 +il89 
-127 -i 11 


+24 +i43 
-320 -i 12 


± 52± 133 


IAM III 


748-i68 


889-i23 


440-i216 


972-i8 


1091-i52 


754-i230 


(errors) 


± 31± i 29 


± 22± i8 


± 7± i 18 


+21, ; 7 


+ 19 +i21 
-45 -HO 


± 22± i 27 



Table 3. Mass and widths or pole positions of the light resonances quoted in the PDG 12 Recall 
that for narrow resonances y/s po [ e ~ M — iF/2 



PDG2004 


p(770) 


K*(892) ± 


cr or / (600) 


/o(980) 


a (980) 


K 


Mass (MeV) 
Width (MeV) 


775.8 ± 0.5 
150.3 ± 1.6 


891.66 ± 0.26 
50.8 ±0.9 


(400-1200)-i (300-500) 
(we list the pole) 


980 ± 10 
40-100 


984.7 ± 1.2 
50-100 


not 
listed 



threshold, which can induce a considerable distortion in the resonance shape, whose 
relation to the pole position could be far from that expected for narrow resonances. 
In addition these states have a large mass and it is likely that their nature should 
be understood from a mixture with heavier states. 



5. Resonance Behavior in the large N c limit 

First of all, let us recall that QCD not only predicts qq states to become bound 
states in the N c —> oo limit, but it tells us exactly what is the large N c dependence 
of their mass and width: their mass should remain constant, whereas their width 
vanishes as l/N c . For instance, the ir,K,rj masses behave as O(l) and fo, their 
decay constant in the chiral limit, as O (yfNr) - 

The N c scaling of the Li parameter d ' ' is listed in Table 1 . Let us then scale 
fo — * fo\/N c /3 and — > Li(p)(N c /?>) for i = 2,3,5,8, keeping the masses, 
2Li — L*2, Li, Lq and L7 constant. We use IAM set II in Table 1 for the parameters 
at N c — 3. In Fig. 2, we show how the poles, represented by a dot, move in the lower 
half of the complex plane as N c changes. On top, we see the p(770) and if* (892) 
vector mesons, whose poles move toward the real axis. That is, the vectors widths 
vanish at large N c , thus becoming bound states. In contrast, in the bottom, we see 
that both the a and k poles move away from the real axis, which means that they 
dissolve in the continuum. 

Not only we can see that vectors become bound states and scalars do not, but 
we can determine what is the precise N c dependence. Before doing that, we should 
nevertheless recall that it is not known at what scale n to apply the large ^V c 
to the Li. Certainly, the logarithmic part in eq.Q) is subleading, but it has been 
pointed ou that it can be rather large for N c = 3, which is our starting point 
for the -ZV C evolution. In addition, the scale dependence is suppressed by 1/N C for 
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Fig. 2. Large N c behavior of poles in the lower half of the second Riemann sheet of the unita- 
rized one-loop meson-meson scattering amplitude from ChPT. For each value of N c the pole is 
represented by a dot, in different meson-meson scattering channel. Note that the a and k behavior 
is opposite to that of well know vector states as the p and K* . 



Li = Li-, L3, L5, Lg, but not for 2L± — L2, £4, L§ and L7. The choice p, or, in other 
words, the choice of "initial values", is a systematic uncertainty typical of the large 
N c approach. Customarily, the uncertainty in the p where the N c scaling applies is 
estimated^ varying p between 0.5 and 1 GeV. We will check that this estimate is 
correct with the vector mesons, firmly established as qq states. In addition, and in 
order to show that the choice of I AM sets in Table 1 is irrelevant for our analysis, 
we change now t o set III. The behavior is exactly that already found in Fig. 2 and 
in previous work^El. 

Let us then look back at vector mesons. In Fig. 3 we show, for increasing N c , the 
modulus of the (I, J) = (1, 1) and (1/2, 1) amplitudes with the Breit-Wigner shape 
of the p and if* (892), respectively. There is always a peak at an almost constant 
position, becoming narrower as N c increases. We also show the evolution of the p 
and K* pole positions, related to their mass and width as ^Js po ie — M — iT/2. We 
have normalized both M and T to their value at iV c = 3 in order to compare with 
the expected qq behavior: Mn /Ms constant and rjv c /T3 ~ 1/N C The agreement 
is remarkable within the gray band that covers the uncertainty ji = 0.5 — 1 GeV 
where to apply the large _/V c scaling. We have checked that outside that band, the 
behavior starts deviating from that of qq states, which confirms that the expected 
scale range where the large _/V c scaling applies is correct. 

In Fig. 4, in contrast, all over the a and k regions the (0,0) and (1/2,0) ampli- 
tudes decrease as N c — ► 00. Their associated poles show a totally different behavior, 
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Fig. 3. Left: Modulus of -kit and nK elastic amplitudes versus y/s for (I, J) = (1, 1), (1/2, 1): 
N c = 3 (thick line), N c = 5 (thin line) and N c = 10 (dotted line), scaled at fj, = 770 MeV. 
Right: p(770) and K*(892) pole positions: ^/s po ; e = M — ir/2 versus JV C . The gray areas cover 
the uncertainty 7V C = 0.5 — 1 GeV. The dotted lines show the expected qq large N c scaling. 
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Fig. 4. Top) Left: Modulus of the (I, J) = (0,0) amplitude, versus ^/s for N c = 3 (thick line), 
N c = 5 (thin line) and N c = 10 (dotted line), scaled at fi = 770 MeV. Center: a mass N c behavior. 
Right: a width N c behavior. Bottom: The same but for the (1/2, 0) amplitude and the ft. 



since their width grows with N c , in conflict with a qq interpretation. (We keep the M, 
r notation, but now as definitions). This is also sugge sted u sing the ChPT leading 
order unitarized amplitudes with a regularization scal c ! ^^ ! 

In order to determine their spectroscopic nature, we abandon momentarily the 
poles and look at the amplitudes in the real axis and the representative Feynman 
diagrams given in in Fig. 5. Imaginary parts are generated from s-channel inter- 
mediate states, as in Fig. 5(a) or 5(c), when they are physically accessible. If there 
was a qq meson, with mass M ~ 0(1) and T ~ 1/N C , we would expect at precisely 
t/I ~ M that Imt ~ O(l) and a peak in the modulus, as it is indeed the case 
of the p and K* . However, we have checked that in the whole a and k regions, 
Imt ~ 0(1/ N%) and Ret ~ 0(1/N C ), which means that from qq states, the a and k 
can only get real contributions from p or K* t-channel exchange, respectively, as in 
Fig. 5(b). The leading s-channel contribution for the k and a comes from Fig. 5(c). 

Then, o ne c an immediately interpret the k as a qqqq (or two meson) state, which 
is predicted^ to unbound and become the meson-meson continuum when N c — > 00. 
However, in the large N c limit, qqqq and glueball exchange count the same, and our 
N c argument alone cannot decide between both structures for the a. Nevertheless, 
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Fig. 5. Representative diagrams contributing to meson-meson scattering and their N c scaling. 



given the fact that glueballs are expected to have masses above 1 GeV, and that the 
K, which has strangeness and thus cannot mix with a glueball, is a natural SU (3) 
partner of the a, a dominant qqqq or two-meson component for the a seems the 
most natural interpretation, although most likely with some glueball mixing. 

Let us finally turn to the other members of the hypothesized scalar nonet: the 
Jo (980) and the ao(980). These resonances are more complicated due to the distor- 
tions caused by the nearby KK threshold, and poles are harder to follow. Still, by 
looking in Fig. 6 at the modulus of the amplitude (0, 0) in the vicinity of the /o(980), 
we see that the characteristic sharp dip of the /o(980) vanishes when N c — > oo, at 
variance with a qq state. For N c > 5 it follows again the 1/iVf scaling compati- 
ble with qqqq states or glueballs. The ao(980) behavior, shown in Fig. 7, is more 
complicated. When we apply the large N c scaling at fi = 0.55 — 1 GeV, its peak 
disappears, suggesting that this is not a qq state, and Imiio follows roughly the 
1/Af? behavior in the whole region. However, as shown in Fig. 5, the peak does not 
vanish at large N c if we take \i = 0.5 GeV. Thus we cannot rule out a possible qq 
nature, or a sizable mixing with qq 7 although it shows up in an extreme corner of 
our uncertainty band. 

Before concluding we want to remark that similar results have been obtained^ 
for all resonances and just for central values of the Li plus a cutoff, using the 
approximated IAM . Furthermore, either the tetraquark structure of these light 
scalar states, or the fact that the lightest qq stat es are above ~1 GeV, has received 
further support from other large N c studieplll 

as well as other kind of analysiJ23. 




Fig. 6. Left: Modulus of a (I, J) = (0, 0) amplitude, versus y/S, for N c = 3 (thick), iV c = 5 (thin), 
N c = 10 (dashed) and N c = 25 (dotted), scaled at ft = 770 MeV. Right: Imt o versus N c . 
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Fig. 7. Left: Modulus of a (I, J) = (1,0) scattering amplitude, versus y/s, for N c = 3 (thick), 
N c = 5 (thin), N c = 10 (dashed) and N c = 25 (dotted), scaled at fj, = 770 MeV. Right: scaled at 
fi = 500 MeV. Center: Imtoo versus N c . The dark gray area covers the uncertainty fj, = 0.55 — 
1 GcV, the light gray area from fi = 0.5 to 0.55 GeV. 



6. Conclusions 

We have reviewed a recent set of workJAI^I^I j n which we show how the unitarized 
one-loop Chiral Perturbation Theory (ChPT) amplitudes generate poles associated 
to the lightest vector and scalar resonances and we study their large N c behavior. 
This amplitudes respect the chiral expansion and are fully renormalized. Indeed 
they provide a remarkable description of two body scattering of pions, kaons and 
etas up to 1.2 GeV using just the 0(p 4 ) ChPT parameters, with values compatible 
with previously existing determinations in the literature. 

We have then studied^ the evolution of the poles, mass and width associated 
to each one of these resonances, through the QCD large iV c scaling inherited by 
the ChPT parameters. We have found that the p and K* vector mesons follow 
remarkably well their expected qq behavior, both qualitatively and quantitatively. 
In contrast, the a, k, /o(980) and ao(980) large N c behavior is in conflict with a 
qq nature (not so conclusively for the ao(980)), and strongly suggests a qqqq or two 
meson main component, maybe with some mixing with glueballs, when possible. 
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